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ABSTRACT 


Many organizations exist within the Department of Defense 
which exhibit the properties of multi-channel, time-dependent 
queueing systems. One such system is the Army's Acute Minor 
Tllness Clinic (AMIC). Models of this system can be developed 
to determine optimum staffing levels and to upgrade the guali- 
ty of service provided. 

This study developed a time-dependent simulation program 
which was applied to a two queue, multi-channel queueing sys- 
tem. In particular, this model was designed to test and 
evaluate the results of an analytical model of the AMIC. The 
results of both of these models are compared to a second 
Simulation program which more closely models the AMIC in order 
to measure the significant differences between it and the two 


other models. 
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I. BACKGROUND 


A. AMIC PROGRAM DESCRIPTION 

An analytical model has been developed utilizing an in- 
homogeneous Markov chain and transition matrices to approxi- 
mate solutions of a multi-channel, time-dependent queueing 
system. The system that was modeled is the Acute Minor I1l- 
ness Clinic (AMIC) which is currently in use by the U. S. Army 
at about twenty of their hospitals. The system integrates 
doctors and specially trained physician assistants called 
Amosists into a dual service process. 

Patients are screened as they arrive by a group of Amo- 
Sists at a reception desk called the Triage. The nature of 
the complaint is determined and the patients are sent toa 
doctor or Amosist for treatment. A patient may request to be 
seen by a doctor in any event. Patients are treated ona 
first come, first served basis as an Amosist or doctor becomes 
available. 

Amosists treat patients according to set procedures de- 
tailed in technical manuals and occasionally on computer pro- 
grams. These procedures, set up by the staff of physicians, 
indicate the actions an Amosist must take in treating a 
patient, and require the Amosist to seek consultations with a 
physician when the problem is beyond his abilities. This pro- 
cedure helps to eliminate errors in diagnosis and treatment. 


When a consultation is indicated, the Amosist locates the 





first available physician and together they treat the physi- 


Cian. Figure 1 depicts the patient flow of this system. 


B. ANALYTICAL MODEL DEVELOPMENT 

Under the direction of Professor Richard W. Butterworth, 
an analytical model of the stochastic process was developed 
with certain simplifying assumptions of the AMIC. Statistics 
were gathered by David L. VanAsdlen and Leonard O. Wahlig, 
and an initial model was implemented. The assumptions of 
this model are listed below. 

First, the arrival stream is an inhomogeneous poisson 
process. Log sheets for 75 days were used to determine 
arrival rate function, A(t), for each day of the week. Upon 
examination, it was determined that the arrival function is 
the same for each day of the week with the exception that the 
expected number of arrivals varies. The function is divided 
into 32 time periods with a mean arrival rate for each seniods 
Figure 2 is a graph of the input arrival function for the 
model. 

Second, the distribution of services is exponential. The 
Amosist system has a distribution of service times for pati- 
ents that do not see a doctor and a distribution of times for 
patients that are referred. The difference in mean service 
times is a result of the amount of time that an Amosist waits 
with the patient for a doctor to become available and to 


treat the patient. The doctors also have two distributions 
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DIAGRAM OF PATIENT FLOW IN THE REAL CLINIC 
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for service, one for referral patients and one for non-refer- 
rals. The mean of the distribution for physician referral 
patients is considerably shorter than that for patients sent 
directly from the Triage. This is due to the nature of the 
consultation. Mean service times are based on 75 observations 
over a three week period. 

Thirdly, in the model the Amosist and the patient do not 
consult with the doctor. Instead, the queues are kept 
separate and independent by adding the percentage of referral 
patients directly to the physician queue and the Amosist 


queue and computing service times as if the consultation took 


place. 
Let: 
a = percentage of patients sent directly to the 
physician queue. 
B = percentage of Amosist patients referred toa 
AOGEOr. ‘ 
Y = percentage of patients that are seen by a doctor. 
Then: 
y=a+t (1 - a)B 
Let: 
E(1) = mean Amosist non-referral service time 
E(2) = mean Amosist referral service time 
E(3) = mean doctor non-referral service time 
E(4) = mean doctor referral service time 
Then: 
(l1-a)8 = percentage of referral patients 
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Mean Amosist encounter time = (1-8) E(1) + (8) E(2) 


Ch) BS), esd) Ba (Ea) 


Meanma@OeCrGrscncOounter tElmMe: mE = ena Seo Se 


In each case the mean encounter time is the weighted sum of 
the average service times. 

Finally, patients may arrive at the clinic prior to the 
time that the clinic opens. This was modeled under the assump- 
tion that the number of arrivals is distributed poisson, for 
each of the two queues. 

Other parameters for the model include the schedules of 
the Amosists and the physicians, and the volume [italics added] 
or total expected number of patients for the particular day 
under study. For detailed studies of the analytical model, 


see References [1] and [5]. 
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it oO uA LON MODE CASE 1 


A. ASSUMPTIONS 

This study developed a simulation program to validate and 
test the significance of the analytical model previously dis- 
cussed. A comparison of resultS gives an indication of the 
correctness of the model's logic and precision by using the 
same assumptions. 

The simulation uses the same parameters as the analytical 
model. It assumes an inhomogeneous Poisson arrival process 
with three types of patients. These types are used to dis- 
tinguish the number of patients arriving at each location. 
Other than for counting and filing in appropriate queues, 
each patient is considered the same by a server. 

Secondly, patients are served uSing a mean service rate 
Midean exponential distribution for each queue. Patients are 
filed in queues upon their arrival and there is no interaction 
between the queues or the servers. Additional variables are 
added to the physician queue for each referral-type patient 
placed in the Amosist queue. 

Thirdly, for purposes of the model, patients who arrive 
prior to the opening of the clinic begin experiencing waiting 
time when the clinic opens. Figure 3 is a flow diagram which 
depicts patient flow for both the analytical model and the 


Simulation model. 
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B. MODEL IMPLEMENTATION 

MiesmModel tStwreteEten in Simscript 22.5, a proprietary 
language supported by Consolidated Analysis Centers, Inc., 
12011 San Vicente Blvd., Los Angeles, California, release 
8D, for implementation on the IBM 360 batch processing system. 
It 1S an event increment program with first in, first out 
routines. The model generates a daily opening and closing 
of the clinic, beginning and ending empty. As such, each day 
1s independent of the others. The simulation is based on a 
particular day of the week for as many iterations as requested. 

Each day the model begins by generating two poisson ran- 
dom deviates. These are the number of patients which arrive 
prior to the opening of the clinic. They are drawn from 
poisson distributions whose means are input parameters to the 
model. A binomial random deviate is then drawn from a dis- 
tribution with the number of trials equal to the number of 
arrivals for the Amosist queue and the probability of success 
equal to g, the percentage of referrals by an Amosist. This 
is the number of Amosist arrivals which are referred to a 
physician. 

Staffing levels change throughout the day as set by input 
parameters to deal with periods of congestion and to allow 
mae StaLt to eat meals. Throughout the day arrivals are gen- 


erated randomly and added to the appropriate queues. 


Ay 





The service times are drawn from exponential distribu- 
tions with one mean for each of the queues as discussed in 


the explanation of the analytical model. 


Mean Amosist service time = (1 -8)E(1) + (E(2))8 


(0) E (3) , (1 -a) 8 (E(4)) 


Mean doctor service time: ee (=a Seoul ane ae 


Cea RRIVAL PROCESS 

The arrival stream is generated uSing the arrival rate 
function supplied to the model. Figure 2 is a graph of the 
function developed from statistics of the Amic. The arrival 
rates for a particular day are computed using the input para- 
meter, volume [italics added]. The volume [italics added] 
is the expected number of patients for a particular day. 
Since part of the daily volume [italics added] arrives prior 
to the opening of the clinic, the arrival rate function does 


not account for all arrivals. Therefore, the arrival rate is 


computed as follows: 


Let : 
V = average number of arrivals per day. 

ave(1l) = mean number of arrivals for Amosist queue prior 
to opening time. 

ave(2) = mean number of arrivals for the physician queue 
prior to opening time. 

A(t) = expected number of arrivals in time [0,t]. 
or A(t) = 


iy (x)dx 
0 
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where 

A(t) = the arrival rate function shown in Fig. 2. 
Scaling the arrival rate function so as to simulate a given 
average volume of arrivals, v, while keeping the same relative 


shape of the arrival function, gives 


A'(t) = v - ave(1l) - ave(2) x A(t) 
ll 
where T is the finite horizon (length of the day). The 


scaled arrival rate function A'(t) is then actually used. 

The actual times of arrival are determined using the cumu- 
iaelvye arrival function. If a set of vat of independent ex- 
ponentially distributed variables are drawn with unit mean, 


Let: 


a! 
S. = .2 
a j= 


ines 


Then S,'s form a homogeneous poisson process. Using the cumu- 
lative arrival function A(t), a relationship exists between 
this homogeneous process and the inhomogeneous process de- 
sired. 

Letting R, = A (S.), the {R,} is one realization of the 
inhomogeneous poisson process and each R, is the time of an 


arrival. Figure 4 is a graph of the cumulative arrival func- 


tion which is computed by the simulation. 


eee SolLATISTICS 
Statistics are accumulated for any specific time of the 


day. Because of the nature of the inhomogeneous process, the 
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average for entire days of system statistics, such as system 
size, do not have a useful meaning. The statistics for each 
parameter are the mean and the standard deviation. 


nN 
we 5° x. SN s=Wcr Oy oe 
: a i <a 


i=l 1= 


The sum and sum of the squares are collected during the simu- 
lation and the means and standard deviations are computed at 
the end. Virtual delay [italics added] is the actual delay 
experienced by each patient during specified time periods. 
The objective is to determine the delay a patient arriving 

at a specific time would experience. In accomplishing this, 
patients who arrive five minutes prior until five minutes af- 
ter a specific time are tagged and their delays are computed 
and averaged. Another statistic, probability that delay is 


[italics 


t 


greater than an input value called "delay criterion,’ 
added} is computed by determining the number of patients in 
the time period that experience delay in excess of that val- 
ue. 

A list of the output statistics follows with definitions: 


Average System Size = Total number of patients in service 
or the queue at a specific time. 


Average Queue Size Total number of patients in the queue 


ata Speci floc tlame. 


Utilization Rate 1 - number of idle servers/number of 


servers on duty. 


(prob that delay=0) 


i-cOeMEeeeSO Sr Vel IS Avallable: @.0 2f 
a server is not available. 
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Virtual Delay = delay experienced by patients at a 
specific time. 


(eazob that delay > = 1.0 if delay is greater than delay. 
delay.criterion) criterion; 0.0 Otherwise. 
E. INPUT TO THE MODEL 

The first part of the simulation includes a routine to 
call input parameters. An example of input parameters is 
appended to this study. Each group of parameters has an op- 
tion word which defines the parameters to be input. The 
routine defaults all parameters to Base Day Monday if input 


is not provided. Option words and their definitions follow. 


Soeron Word Definition Parameter 
ENCOUNTER Mean service l-dimensional array > = 0 
time Encounter (1)=Amosist nonreferral 


mean service time 


Encounter (2)=Amosist referral 
mean service time 


Encounter (3)=doctor nonreferral 
mean. censyvice time 


Encounter (4)=doctor referral 
mean service time 


Default values: 


Encounter (1) = 17.96 
(2) = 23.83 
(3)0= 13.55 
(4) = 4.34 
VOLUME. Average num- Integer > 0 
meee ana cal Default value: 
a Y Volume = 156 
OPENING. 24-hour clock real > 0 
TIME time that Default value: 
clinic opens Opening.time = 06800. 


22 





CLOSING. 
TIME 


PERCENTAGE 


REFERRALS 


ITERATIONS 


PRELOADING 


AMOSIST.SCH 


PHYSICIANS 


24-hour clock 
time that 
clinic closes 


Percentage of 
patients tri- 
aged directly 
EO ca saOC LOE 


Percentage of 
patients re- 
ferred by an 
Amosist 


Number of days 
of simulation 


Average number 
of patients 
erigied Vag. (Pi OT 
to opening 


Schedule for 
Amosist staff 


Schedule for 
Physician staff 
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real > 0 
Default value: 


Closing.time = 2400 
real > OQ 

Default value: 
alpha = .55 

Bea 240 

Default value: 

Beta = .3l 


integer > 0 
Default value: 
Last.day = l 


1-dimensional array 

Ave(1l) = number of Amosist 
arrivals 

Ave(2) = number of doctor 
arrivals 

Default Values: 

Ave.arrival (1) 

Ave.arrival (2) 


> ~] 


integer > 0 

No.duty.chgs = number of staff 
changes 

Hour(i) = time shift starts 

No.amosists(i) = no of staff 
fOr yeEhat sshitt 

Hour(1l) must = opening.time 

Default values: 

NOsGiney -chas7— 5 


Hour(1l) = 0800 No.amosists(1)=5 
(2) = 1100 (2)=3 
(ee — e300 (3)=5 
(4) = 1530 (4)=7 
(S)> = 1630 (5)=2 


integer > 0 
No.rosterchgs = number of staff 


changes 
Clock(i) = time shift starts 
No;-C@OecEOrsia) = number of staft 


PObLethac SNILe 
Clock(1) must = opening.time 
Default values: 
Clock(1)=0800 No.doctors (1)=2 
(2)=1100 (2)=1 





ARRIVAL. FCN 


DELAY . VALUE 


Arrival function 


Time in minutes 

to compute prob 

delay is greater 
than delay 

SE heciaron 


24 


(2) 
= 
(2) = 


integer 
pairs = 


Time (i) 


Default 
Pairs = 
time (1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
(11) 
(ee?) 
(il3) 
(14) 
(15) 
(16) 
(ky) 
(18) 
(19) 
(20) 
(2) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 
(29) 
(SO) 
(315) 
(eZ) 


real > 0 
Default value: 
Delay.criterion 


1200 
is 630) 
1630 


e0) 


(22 
(ao 
ia 


no of periods de- 
fined 
2,1-dimensional arrays 
= end of period 
Customer(i) = 
arrival 
alternate time and customer 


rate of 


values: 
82 
= 810 Customer (1)=22.5 
= 820 (2) =23.5 
= 840 (3)=24.5 
eo) (4)=23.5 
= 900 (Sy 22.5 
= 910 (6)=21.5 
= 920 (7) =20 25 
= 930 (8)=19.6 
= 945 (9)=18.6 
=1000 (10) =1.7 26 
= INOS (11)=16.6 
— 0G 102.) wea 
=1045 (13) =13.35 
=1100 (14)=12.05 
=1115 (15) =00.95 
=1130 (16)=9.9 
see Ch} =8.9 
=1200 (18)=7.9 
— ees (19)=7.0 
1245 (20)=6.2 
=1300 (21)=7.2 
=1330 (22)=13.6 
=1400 (25) —l2e5 
=1430 (7725) = eel 
=1500 (25)=9 285 
=o O (26)=8.8 
=1630 (27) =7~.55 
=e 0 (QS =6.5 
=1930 (29)=5.8 
= eiG (30)=4.4 
=2330 (31h —2 38 
=2400 (32)=.01 
= 15.0 





END.OF.DATA Bnd Of Input. Nothing should follow this 
Cana. 

The input routine is designed to be as user oriented as 
possible. All input is in free format with a blank space be- 
mage tne deliminater for each input word. The input option 
cards and their associated values may be read in whatever 


order desired as long as the end.of.data card is last. 


F. VERIFICATION OF THE PROGRAM 

In the early stages of formulation extensive output was 
used to follow, exactly, the patient flow through the model. 
In its final form the model gives a daily output of the nunm- 
ber of patients treated in each system by type, the number of 
arrivals prior to opening time for each queue, and the time 
the clinic closed on that day. Provision has been made, if 
the congestion continues throughout the day, to count the 
number of patients that are waiting for service when the 
clinic closes. In the actual clinic these people would be 
Sent to the emergency room for treatment. 

The number of referred patients served in both queues was 
always equal and the sum of the referral patients and the 
non-referral patients from both queues always equalled the 
total number of arrivals for the day. 

The means of the poisson deviates and the binomial devi- 
ate from the early arrivals were determined. For patients 
seen only by a doctor the mean was 7.2 with an input mean of 


7 for a 200 day simulation. For the amosist queue the mean 
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was 3.92 with a standard deviation of 1.95 and an input para- 
meter of 4. The expected number of referral patients was 2.17 
and the actual average was 2.24 with a standard deviation of 
1.47. There is no significant difference in the means. 

The total number of arrivals for each queue was also 
summed by type over the simulation to compare to the expected 
values, to insure the correct percentage of patients was 
directed to each queue. 


Ave. num. of doctor's patients = 84.16; expected 
num = 85.8 


Ave. num. of amosist's patients = 48.94; expected 
num = 48.44 


Ave. num of referral patients = 22.36; expected 
num = 21.76 


A test of hypothesis was made on the arrival distribution 
to determine if the number of arrivals for each period was 
from the same distribution. The statistic is Q and is dis- 


tributed vane) ‘ 


8) 


2 (xX; - nr.) 
arr), 


: ig aie 
le 


1 


where rs is the mean number of arrivals in the ith period and 
x. is the observed number of arrivals over all iterations. 
At the .025 level of Significance the ne is 47.0 with 31 de- 
grees of freedom. The test statistic was 37.64, and the hypo- 
thesis was accepted. 

The mean service times were computed using a language 


generated routine, and compared to the expected values. Since 
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the program uses internal exponential generating routines, it 
was only necessary to check the means. 


Amosist mean service time = 19.583; expected value 19.78 


Doctor mean service time 11.559; expected value 11.687 


The results of the model are tabulated ina later section 
of this study. Source listings and flow charts of this pro- 


gram are appended to the study. 
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IIt. CASE 2 SIMULATION MODEL 


A. ASSUMPTIONS AND CHANGES TO THE MODEL 

After developing a model of the system using the assump- 
tions of the analytical model, the next step in the study was 
to produce a model whose assumptions more closely duplicated 
those of the clinic. The change required was in the handling 
of referral patients. In this model as in the actual system, 
all patients going to the Amosist have one distribution of 
service times. In the previous model, the extra time for a 
referral patient was time used waiting to see a physician and 
consulting with him. The Amosist is held out of service in 
this model until the consultation is complete and then the 
Amosist can continue treating other patients. In the physi- 
cian queue, referral patients preempt non-referral patients 
and have a separate distribution of service times from the 
patients who are sent directly from the Triage to a doctor. 
All other assumptions of the model remain the same as the 
Case 1 model. Figure 5 is a flow diagram for patients in 


this model. 


Bees VERIFICATION 

The input and arrival process for this model was unchanged 
and, therefore, was not retested. The service times averaged 
18.119 minutes for the Amosist queue and 11.419 minutes for 


the physician queue. The expected values were 17.96 minutes 


PAS) 








DIAGRAM OF PATIENT FLOW_IN THE 
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Amosist service channels during a consoltation. 


Figure 5. 
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for the Amosist queue and 11.687 minutes for the physician 
queue. A source listing for this model is appended to the 


study. 
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Peo Una One RESULTS 


A. COMPARISON OF ANALYTICAL AND CASE 1 MODELS 

In order to test the results of the two simulations against 
the analytical model of the AMIC, the means and standard de- 
viations of each of the output parameters were tabulated. The 
results of the analytical model are found in Table I. Table 
II is the results of the Case 1 simulation model. 

Since the simulations and the anlytical model have 
separate distributions and alternate methods of obtaining 
values which are independent, a test can be made that the 
average values are now significantly different. 


Ho: X = A where X 
A 


Simulation statistic 
analytical deterministic value 


For purposes of the test the output parameter of the analyti- 


cal model is assumed to be the true parameter of the system. 


mec > 
N = number of samples in the distribution (2)) 
X = mean value of simulation parameter 
A = mean value of analytical model 
: = standard deviation of X 
T = test statistic 
Then: 


ey SS Te yy 
S 
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RESULTS OF THE CASE 1 SIMULATION MODEL 
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RESULTS OF THE CASE 1 SIMULATION MODEL 
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Since there are 200 iterations in the simulation, the T sta- 
jmMatrtc approaches its limiting distribution which is normal 
(0,1). Tables III and IV are the summary of this comparison. 
All "Undef" entries indicate a standard deviation of 0.0 
indicating that all values of X were equal at that time. 

The two models were compared for the parameters of base 
day Monday. A critical region for accepting or rejecting the 
hypothesis was not chosen since entries are not independent 
of time. Although each day's statistics are independent, the 
Statistics taken at different times of the day are serially 
correlated, and cannot be pooled. Therefore, subjective 
judgment based on the normal statistics computed was used to 
judge the fit of the two models. 

Using the above criteria, the analytical model is a very 
good fit in the Amosist system. However, when congestion 
overloads the queues, the analytical model does not conserve, 
pamper che probability used in computing the output para- 
meters. As a result, the analytical model underestimates 
the number of patients in the system. Because the physician 
system in this particular set of parameters was highly con- 
gested, the physician system of the analytical model indi- 
cated fewer patients being serviced. The values of both 


models were very close after the congestion diminished. 
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TABLE IIf,. COMPARISON OF THE ANALYTICAL AND SIMULATION MODELS 


Megas octal ISliCs FOR THE CASE J AMOSIST SYSTEM 


System size Queue size Utilization Prob(Del=0) Average Delay Time 
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TABLE IV.  COMPARTSON OF THE ANATYTICAL AND SIMULATION MODEUS. 


MORMAL, StALTISTICS FOR THE CASE 1 PHYSICIAN SYSTEM 


System size Queue size Utilization Prob(Del=0) Average Delay Time 
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B. COMPARISON OF THE CASE 2 SIMULATION AND ANALYTICAL MODELS 

The same procedure outlined for Case 1 was used in com- 
paring the results of the Case 2 model. Table V is the re- 
sults of the Case 2 model and Tables VI and VII are the sum- 
mary of the comparison tests for the models. 

In this case other factors caused differences between the 
two models. In the Case 2 model, Amosists spend more time 
than expected waiting for consultations. Therefore, the two 
models have different distributions for referral service. 

The result was that the Amosist queue was longer than the 
Original model and the analytical model results were not as 
close as previously seen. In the other system, the physician 
queue was smaller than in the Case 1 model. This was the re- 
sult of the phenomenon of shortest processing time. Since all 
of the patients with the shorter processing distribution were 
served first, the size of the queues are expected to be short- 
er. The average time in the system for a patient does not 
change, but because of the preemption fewer people are wait- 


ing for service at any one time. 


C. COMPARISON OF CASE 1 AND CASE 2 MODELS 

As discussed in the previous section, the average service 
time for an Amosist referral patient was longer than expected 
in the Case 1 model or the analytical model. The physician 
system had shorter queues and fewer patients in the system 
at any one time due to the fact that the distribution of 


service times for referral patients was much shorter than the 


39 





-_— 











RESULTS OF THE CASE 2 SIMULATION MODEL 
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RESULTS OF THE CASE 2 SIMULATION MODEL 
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Peete Vig COMPARISON OF THE ANALYTICAL AND SIMULATION MODELS 
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TABLE VIT. COMPARTSON OF THE ANALYTICAL AND SIMULATION MODELS 


Perv SAP lot tGS FOR THE CASE 2 PHYSICIAN SYSTEM 


System size Queue size Utilization Prob(Del=0) Average Delay Time 
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ester i bution for other patients. The result of servicing 
these shorter processing times first, causes the shorter 
queues. However, the output statistics of the two models were 
quite close, indicating that the Case 1 model is a good model 
of the real system. The fact that the output deterministic 
values for the analytical model were so close to the other 

two models is an indication that this analytical model is a 


good approach to solving the problem. 
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V. CASE 3 SIMULATION MODEL 


foe HOSUMPTIONS 

The weakest of the assumptions in all of the models was 
the assumption of exponential service times distributions. A 
third simulation was conducted to determine the sensitivity 
of the model to the service distributions. A degenerate dis- 
tribution was used for these service times. The constant 
value was equal to the means of the exponential service dis- 
tributions used in the other models. All other assumptions 
remained the same as in the Case 1 model. Table VIII is a 


Summary of the results of the simulation. 


Dae skcoULTS 

Based on this one variation to the base case, the exponen- 
tial assumption is apparently robust. Due to the randomness 
Meecent in simulations the queves were larger than the previ- 
ous models in the earlier part of the day and smaller in the 
later part of the day. There is not sufficient difference to 
indicate that the exponential distribution assumption is in- 
valid. The means of the service times is critical, however, 


as waS indicated in the previous section. 
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RESULTS OF THE CASE 3 SIMULATION MODEL 
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VI. CONCLUSIONS 


A. COMPARISON OF THE ANALYTICAL MODEL TO THE SIMULATION MODELS 

The purpose of this study was to devise a simulation for 
multi-purpose queueing systems which could be applied to many 
varied systems. Second, the purpose was to determine the pre- 
cision and correctness of an analytical model which had been 
applied to such a system. A first case study, duplicating the 
assumptions of the model indicated that, although some weak- 
nesses existed in the analytical model, the application was a 
good approximation of the AMIC Clinic. The model has a finite 
queue and, therefore, underestimates the system at very high 
levels of congestion. Also, in the real system, the preemp- 
tion of patients by patients whose average service times 1s 
so much less causes smaller queues than a first in, first out 
system. 

The second simulation program revealed the difference be- 
tween the real system and the analytical model in that the 
real world referral patients preempt patients waiting to see 
physicians, and the Amosists waiting with the patient cannot 
service other patients. The distribution of service times 
for these Amosist referral patients has a larger mean than 
Originally believed, causing Amosist statistics for the 
second program to be larger. This program also indicated 


smaller physician statistics due to the preemption. 
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Eee SENSITIVITY OF THE MODEL 

A third simulation model indicated that a service distri- 
bution assumption of exponentially distributed service times 
was apparently robust. However, the system is very sensitive 
to the means of the service distributions. In all the analy- 
tical model appears to be a good fit, and useful in providing 
the user with information to upgrade service and optimize 


staffing levels. 
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VI. GENERAL NATURE OF THE SIMULATION MODEL 


The simulation programs were developed in as general a 

Way aS possible to make them applicable to many varied sys- 
tems. By manipulating the input parameters, any queueing sys- 
tem with up to two queues and some limitations on interaction 
is capable of being modeled by these programs. For example, 

a system with only one queue 1s possible by setting the per- 
centage input parameter to 1.0, thus sending all arrivals to 
one queue. Several assumptions remain. Services are distri- 
buted exponentially, and the arrivals make up a poisson pro- 
cess. The preloading or allowing of arrivals prior to open- 
ing may be surpressed by setting the average arrival prior to 
opening ae 0.0. The arrival process may be homogeneous or in- 
homogeneous depending on the functional values which are in- 
put. Supplying one arrival rate for the entire day produces 

a homogeneous poisson process, with a mean arrival rate satel 


to the input parameter. 
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